Odd momentum pairing and superconductivity in vertical graphene heterostructures 
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Vertical graphene heterostructures made up of graphene layers separated by BN spacers allow for 
novel ways of tuning the interactions between electrons. We study the possibility of electron pairing 
mediated by modified repulsive interactions. Long range intravalley and short range intervalley 
interactions give rise to different anisotropic phases. We show that a superconducting state with 
gaps of opposite signs in different valleys, an odd momentum pairing state, can exist above carrier 
densities of 5 — f x 10 13 cm -2 . The dependence of the transition temperature on the different 
parameters of the devices is studied in detail. 

PACS numbers: 



INTRODUCTION 

Graphene has a number of electronic properties which 
can be tuned [TJ [5]. Changes in the carrier concen- 
tration, nature of the substrate, and other properties 
of the environment can modify the electron-electron 
interaction^ 0]. Graphene structures, made by com- 
bining graphene layers separated by hexagonal boron 
nitride [5TI10] show a number of interesting features, re- 
lated to the interactions between electrons in different 
layers. These devices are highly tunable, as the carrier 
concentration in the graphene layers, and the distance 
between them can be varied independently. 

We study here the possibility of superconductivity in 
these devices, arising from a suitably modulated electron- 




FIG. f : Sketch of the heterostructure considered in the text. 
Two graphene mono- or bilayers are embedded into three 
dielectric media. The three dielectrics need not have the 
same dielectric constant. The device is characterized by the 
carrier densities in the graphene layers, p t and pt, the dis- 
tance between layers, d, and the three dielectric constants, 
ei = et, £2 = e m and e 3 = e b . 



electron interaction. Superconductivity in metallic mul- 
tilayers due to electron-electron interactions was pro- 
posed long agojll], as these systems exhibit an acous- 
tic plasmon, which can mediate an attractive interac- 
tion between the electrons. The model has been gen- 
eralized to graphene near a metal surface [T^]. Alter- 
natively, a repulsive electron-electron interaction, if its 
spatial dependence is suitably modulated, can lead to 
superonductivitv|13j . This possibility has been exten- 
sively studied in many materials (see, for instance 1 141- 
These studies include carbon nanotubes[2"4"]. 
graphene and graphite at certain dopings or stacking 
arrangements pi, 25-33 , and combinations of two dimen- 
sional electron gases [T51|201 12"3] (2DEG). 

Experimentally, a number of graphene based com- 
pounds exhibit superconductivity. The graphite interca- 
lation compounds [3"4"H3"5] show superconductivity up to 
temperatures T c s» 11K, and carrier concentrations in 
the graphene layers of n ~ 1 — 4 x 10 14 cm~ 2 . A number 
of doped fullerenes have a superconducting phase [39| [40] 
with the highest critical temperature T c ~ 38K. The car- 
rier concentration at which superconductivity is found, 
normalized to the number of carbon atoms, is simi- 
lar to the concentrations in the graphite intercalation 
compounds. The origin of superconductivity in doped 
fullerenes at high pressures is not completely understood, 
and it is likely that electron-electron interactions play a 
rolegl]. 

The carrier density and interlayer distance in graphene 
heterostructures determine the features of the electron- 
electron interaction, which can be changed over a wide 
range. This opens the possibility of superconducting 
phases at not too large carrier densities and without the 
requirement of special features in the density of states. 
The analysis of superconducting instabilities is the main 
goal of this work. The next section discusses in detail the 
effective interaction. Then, we analyze possible super- 
conducting instabilities using the Kohn-Luttinger frame- 
work, and we also characterize the most likely supercon- 
ducting phase. We find that an odd momentum paring 
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state, with gaps of opposite signs in the different cones, 
is favored by intervalley scattering at large densities and 
could be the leading instability. We present next the de- 
pendence of the transition temperature on the main pa- 
rameters of the model, the nature of the superconducting 
phase, and also the relation to the formation of an ex- 
citonic condensate in an artificial bilayer with electrons 
and holes. 



EFFECTIVE INTERACTION IN GRAPHENE 
HETEROSTRUCTURES 

We consider two graphene layers, each of which can 
be single layer graphene or bilayer graphene with carrier 
concentration ro t and rib, where t and b denote the top and 
bottom layers. They are separated at a distance d. The 
top, middle, and bottom dielectrics have static dielectric 
constants £1,62 and £3, as shown in Fig. [T] The param- 



eters used are: nearest neighbor hopping 70 ~ 2.7 eV, 
distance between nearest neighbor carbon atoms a ~ 1.4 
A, Fermi velocity Hvf = (37oa)/2, fine structure constant 
of graphene a = e 2 /(hv F ) « 2.5, and hopping between 
nearest neighbors in different layers 71 « 0.4 eV. 

We assume that Cooper pairs are formed by elec- 
trons which occupy states related by time inversion 
symmetry |42j. We consider processes in which the 
Cooper pairs are scattered within one valley, and between 
valleys, as shown in Fig. [4] We study first the effect of 
the scattering of Cooper pairs within one valley, shown in 
Fig. [4] a) . We approximate the long wavelength electron- 
electron interaction between the layers by its RPA value, 
see Fig. [5} This approximation is justified when the num- 
ber of fermion species, Nf, is large. In our case, Nf = 4. 



The effective interaction for small momenta, 
in the top layer is given by 
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where v c {q) = 2ire 2 /q is the Coulomb potential, and 
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These expressions [43] interpolate the effective dielec- 
tric constant between lim g( j->o £« (q) = (£1 + £3)/% to 
liniqd^oo ett(q) = (£1 + ea)/2. The intra valley static 
polarizations of single layer and bilayer graphene are 
known [441446J . For bilayer graphene we use the expres- 
sion which includes contribution from the four bands [46], 
The electronic excitations of the system, and their cou- 



pling to the electrons in each layer, can be described by 
Imvu(w,q) and Imvb b (uj, q). These quantities are shown 
in Fig. [2] and Fig. [3] for different values of the screening 
by the dielectric layers. The calculations include a finite 
broadening, in order to display the plasmon resonance. 
The charged plasmon has a large fraction of the spec- 
tral strength, and its dispersion depends on the screening 
constant. The acoustic plasmon has a dispersion 
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For the parameters used in Fig. [2] and Fig. [3] the plas- 
mon velocity is very close to the Fermi velocity, and its 
spectral weight in the effective interaction is small. 
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KOHN-LUTTINGER SUPERCONDUCTIVITY IN 
GRAPHENE HETEROSTRUCTURES 

Pairing interaction 

Superconductivity due to electron-electron interac- 
tions is proportional to the Fermi energy, tp, times 
a factor which depends exponentially on the inverse 



of the coupling, A, T c ps X.ltFe -1 ^"-' . In materi- 
als where superconductivity is induced by phonons, the 
Fermi energy is replaced by the Debye frequency, T c pa 
l.lhd£,e~ 1 / X ' l -p h . Because of the difference between the 
Fermi and Debye energies in a given material, electron- 
electron interactions have been considered a route to 
high critical temperatures, provided that the couplings 
are similar. The effective electron-electron interaction 



FIG. 2: Imaginary part of the effective potential for an heteroestructure with p t = 4 x 10 13 cm 2 ,p& = 10 13 cm 2 ,ei = 62 
£3 = 4, and d = 10 nm. Left: top layer. Right: bottom layer. 
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FIG. 3: As in Fig. [2] but with e< = 10. Left: top layer. Right: bottom layer. 



in an isolated 2DEG, however, is only weakly modu- 
lated as function of momentum q in the relevant range 
< |q| < 2kp, where is the Fermi momentum|47j. and 
the Kohn-Luttinger instability, although present, leads to 
low critical temperatures. 

In a two dimensional metallic system with an isotropic 
Fermi surface, the coupling constant A„ for a supercon- 
ducting phase with an order parameter A„(0) oc cos(n#), 
is given by 

2p(e F ) 
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where p{cf) is the density of states per spin at the Fermi 
level, and V^(q) is the effective electron-electron interac- 
tion. For superconductivity to occur, A„ must be neg- 
ative. The critical temperature is given, approximately, 
by 
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The calculation of A„ is simplified if the function V(q) 
has singularities or a non analytical dependence [13] on 
q = |q|. For instance, when 
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we obtain 
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Hence, superconductivity at large angular momenta is 
possible when V'(0) > 0. In general, superconductivity 
is favored if V(\q\ = 2k F ) > V(0). 



Pairing due to intravalley scattering 

We assume that the distance between the chemical po- 
tential and the edge of the gap induced at the Dirac en- 
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FIG. 4: Sketch of the scattering processes of Cooper pairs 
considered in the text: a) Scattering within a given valley, b) 
Scattering between valleys. 
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FIG. 5: Diagram included in the calculation of the effective 
interaction, V(q). Red and blue bubbles stand for the polariz- 
abilities in the two layers in the heterostructure. The dashed 
lines describe the Coulomb interaction, which can be intra- 
or interlayer. 



ergy is larger than the gap itself, so that the pairing in- 
teraction does not couple the conduction and the valence 
bands. The densities of states per spin are 
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FIG. 6: Effective interaction, A(6) = 

{2/n)p(e t F )V[2k t F sm(9/2)] as function of 9 for an het- 
erostructure made of two single layer graphene sheets with 
densities p t = p b = 10 13 cm~ 2 and d = 10, 40, 70, •■• 280 
A(from bottom to top). 
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FIG. 7: Values of A„ as function of k F d for the same param- 
eters as in Fig. [6] 



The chiral nature of the wavcfunctions in graphene 
changes the expression in eq. [3] into 
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A constant intervalley potential, V(K — K') does not 
contribute to pairing in the situation considered here, 
where the gap is modulated within each valley. Higher 
order terms in an expansion in momenta around K — K' 
fix the relative phase of the gaps in the two valleys. 

As mentioned earlier, pairing is favored if the interac- 
tion potential is such that V(2kp) > V(0). A finite value 
of V(0) is ensured by the metallic screening induced by 
the graphene layers. The value of V(0) is further reduced 
if the dielectric constant of one of the insulating layers is 
much larger than the other two, i. e. e 3 ^> e 1; e 2 , a situa- 
tion which seems feasible in current experiments |481 . 
Then, the screening in one layer changes from (e 3 +ei)/2 
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FIG. 8: Effective interaction, A{9) for a graphene heterostruc- 
ture made of single layer graphene and bilayer graphene with 
densities p s i g = pu g = 10 13 cm -2 and d — 10 A. 




to (e 2 + £i)/2, as discussed before. For two graphene 
monolayers, the values of the dimcnsionless constants 
7„ are determined by a — e 2 /vF,k F d,k F d,ei,€2 and 
£3. We present in Fig. [6] calculations of the function 
A(9) = 2/^p t (e F )V[2k F sin(0/2)] cos 2 (9/2) for p t = p b = 
10 13 cm _2 ,e 1 = e 2 = e 3 = 4, and different values of d. 
The value of A(0) raises for small 0, with a slope which 
increases with k F d. The orthogonality between wave- 
functions of opposite momentum leads to A(n) = 0. 

The coupling is enhanced using insulating layers with 
different dielectric constants. For ei = = 1 and £3 = 6, 
and the densities and distances shown in Fig. [6] we find 
coupling constants of up to A3 ~ —0.055. The critical 
temperatures are of order 
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These temperatures are too low to be observed experi- 
mentally. 

These estimates do not change much for an het- 
erostructure made of single layer and bilayer graphene. 
Fig. [8] shows results for a device with similar parameters 
as in Fig. [6j 



FIG. 9: Sketch of the gap with alternating signs in different 
valleys induced by intervalley scattering. The odd momen- 
tum pairing state, A(k) = — A(— k), breaks time reversal and 
inversion symmetry. 



The value of the on site repulsion term in carbon 
7r orbitals, U, has been studied for a long time, by 
fitting model hamiltonians to the spectra of organic 
molecules [501452] . The value of U is also a basic input in 
the analysis of excitonic peaks in carbon nanotubes |531 
151] . Recent estimates (551 US] suggest that a reason- 
able value of U is U « 10 eV. This value is, approx- 
imately, U « e 2 /a, or, alternatively, Uft « 27re 2 /K, 
where fl = 3V3a 2 /2 is the area of the unit cell, and 

K = 4tt/ (3%/3a). The value U « 10 eV is also consis- 
tent with the value U ~ 11 eV used to fit the position of 
excitonic peaks in carbon nanotubes 53, 54J. 

We neglect the dependence of U on momentum away 
from Q = K K'. Then, pairing is only possible if A^(#) 
and Ak'(0) are independent of 0, and Ak = —Ak<- The 
coupling constant is 
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Pairing due to intervalley scattering 

We consider now the scattering of Cooper pairs be- 
tween valleys, see Fig. [4] b). We take V(K - K') w UCl, 
where U is the Hubbard repulsion between electrons on 
the same carbon atom, and f2 = 3\/3a 2 /2 is the area of 
the unit cell, and a is the distance between carbon atoms. 

The next leading coupling, V , the repulsion Intervalley pairing is favored by the term 

between electrons in atoms which are nearest neighbors, 

does not contribute to the scattering process depicted in jjfi 3y/3Ukpa 2 y/SUkp 

Fig. [i] b). Repulsive interactions between electrons in p( € f)~^~ — ^ Vp = 2^70 

atoms at larger distances can be approximately included 
as an increase in the value of the Hubbard repulsion. The 
combination of weak intravalley and strong intervalley re- 
pulsion leads to a gap with alternating signs in different 
valleys, as sketched in Fig. [9] We assume that the gap 
has a constant value within each valley. 
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so that intervalley pairing is more likely at high densities. 

We can obtain an order of magnitude estimate of the 
strength of the screened intravalley repulsion. In the 
limit kpd — > 0, and assuming a single dielectric constant, 
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FIG. 10: Dependence of the parameter Ap in eq. |16| on distance 
for an heteroj unction with t\ — £2 = £3 = 4, nt = 10 14 cm -2 
and n b = 2 x 10 14 cm~ 2 . 




FIG. 11: Dependence of the critical temperature on carrier 
density for an heterostructure made of two single layers with 
the same concentration, separated by a distance d = 10 A, 
and £i = £2 = £3 = 4. 



60) we obtain 
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where N F = 2 is the number of layers, and Nf = 4 is the 
number of fermion species. When comparing the first 
and second term on the r. h. s. in eq. |11| one needs to 
take into account that the absence of backward scatter- 
ing, described by the factor cos 2 (0/2) in the integrand of 
the intravalley contribution reduces by a factor of about 
1/2 that term. From eqs. [l2| and [13] one can obtain an 
estimate of the range of parameters for which supercon- 
ductivity is possible 
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Approximating U ~ e 2 /a, we obtain a minimum density 
above which superconductivity will appear 
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8 x 10 13 cm 2 . The value of the intervalley pairing term 
for n = 1.6 X 10 14 cm- 2 , is p(e F )Un 0.32. 

A more precise analysis of the effects of intravalley in- 
teractions is described by the parameter 
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The dependence of A and Aq on distance for a het- 
erostructure made from two graphene single layers, and 
ei = e 2 = e 3 = 4,n t = 10 14 cm -2 and n b = 2 x 10 14 cm~ 2 
is shown in Fig. [10] The critical temperatures of the two 
layers are T\ w 0.6K, T b c w 100K. 

The value of T c is highly dependent on the input pa- 
rameters, and is enhanced as the interlayer distance is 
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FIG. 12: Dependence of T c on interlayer distance and dielec- 
tric constant. Top: Dependence on distance, for U = 10 eV 
and £i = £2 = £3 = 4. Bottom: Dependence on e ly for (7 = 10 
eV, and d = 10 A. 



reduced, or else in the presence of a substrate with a high 
dielectric constant, such as HfC>2, see Figs. [Tl] and [12} 

The mean field analysis reported here does not include 
phase fluctuations of the order parameter. In a two di- 
mensional superconductor these fluctuations lower the 
critical temperature, and turn the phase transition into 
a Kosterlitz Thouless transition. The change in critical 
temperature can be accounted for as a correction by nu- 
merical factor of order unity of the mean field value [57], 
The microscopic parameters which define the model have 
a weak influence on this factor. 
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FIG. 13: Dependence of T c on the electron phonon coupling, 
g p h- The remaining parameters are U — 10 eV, d = 10 A, 
ei = £2 = £3 = 4. 



and has opposite signs in the two valleys. Elastic scatter- 
ing within each valley does not change this phase, while 
intervalley scattering is pair breaking. Defects which in- 
duce intervalley scattering induce Andreev states inside 
the superconducting gap. As an example, we consider a 
local perturbation, described by a shift of the onsite en- 
ergy of a single it orbital, which we denote 5e. The local 
density of states in the presence of the perturbation is 



where 
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Electron-phonon coupling 

The leading contribution to the electron-phonon cou- 
pling in graphene is due to optical phonons at the T 
and K points of the Brillouin zone. These phonons 
can induce intra- and intervalley scattering, and opti- 
cal phonons modify strongly the interatomic distances, 
and the nearest neighbor hopping elements. The electron 
phonon coupling induces an attractive interaction at all 
wavevectors, and favors the existence of a superconduct- 
ing gap with the same sign in the two valleys [251 EH- 
Expanding around the Dirac points, the coupling can 
be parametrized by effective intra- and intervalley po- 
tentials, 9p hl gp h , which can be estimated from LDA 
calculations [58], or extracted from experiments |59[ 160] , 
A tight binding estimate of this coupling 60 62J is 
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where Mq is the mass of the carbon atom, u)q,i is the 
frequency of the phonon, and (3 = 91og(7 )/91og(a) ps 
2 — 3. This estimate gives a dimensionless coupling, for 
/3 = 2, huj n = 0.2 eV, which is doubly degenerate, and 
toi = 0.17 eV, 
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The value of this coupling is in reasonable agreement with 
LDA calculations GEO, where the numerical factor is 5.5. 

2<&+ 



The total coupling obtained from eq. 
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gh h ~ 2.3 eV A 2 . This coupling reduces the effect of the 
Coulomb interaction. The resulting effect on the value of 
T c is shown in Fig. |13| 



Andreev states and pair breaking effects in the 
superconducting phase 

The superconducting phase studied in the last sub- 
section shows a gap which is constant within each valley 
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where p(e F ) = 2kp /{irvp) x fl is the density of states at 
the Fermi energy and f2 is the area of the unit cell. Using 
Ak = —Ak>, we find that the Green's function in eq. 
has poles at 
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A vacancy can be described as the limit 5e — > oo. In 
this limit, we find two midgap states, at u = 0. Im- 
purities which induce intervalley scattering induce pair 
breaking and reduce the gap. Superconductivity will be 
completely suppressed for A w vf/£k-K', where Ik-K' 
is the elastic mean free path associated to intervalley 
scattering. 

In a similar way as in other multiple band materials, 
Andreev states can be induced at surfaces [53]. Alterna- 
tively, surface states can be seen as arising from topo- 
logical features of the bulk band. In order to describe 
effects due to the breaking of translational symmetry, it 
is convenient to write the Bogoliubov-de Gennes hamilto- 
nian in a real space basis. A superconducting gap which 
changes sign between valleys in the honeycomb lattice 
can be described by a tight binding model with purely 
imaginary hopping terms which mix electron and hole 
states at sites the same sublattice, following the same 
method used in the formulation of the Haldane model [5^]. 
The tight binding model contains four entries per lattice 
site. In the absence of magnetic impurities, the model 
can be split into two equivalent systems. In the presence 
of electron-hole symmetry, with the chemical potential at 
the Dirac point, each of these models can be further split 
into two Haldane models with opposite Chern numbers. 
Thus, the odd pairing considered here, in the presence 




FIG. 14: Subbands of a graphene ribbon with a supercon- 
ducting gap which has opposite signs in the two valleys. In 
order to illustrate better the main features, the Fermi energy 
is e_F = 0.5 eV, and the gap is A = 0.2 eV. Top: armchair 
edge. Bottom: zigzag edge. 



of full electron-hole symmetry can be seen as the sum 
of independent topological superconductors. For finite 
chemical potentials, however, this reduction cannot be 
done, and the superconducting phase has no topologi- 
cally protected Andreev states at any edge. 

In the absence of general topological arguments, An- 
dreev states can be expected only at surfaces where the 
valleys are hybridized. Fig. [14] show the subbands of su- 
perconducting nanoribbons with different orientations. 
In agreement with the previous analysis, subgap states 
exist for armchair nanoribbons, while they are absent in 
zigzag nanoribbons. Note that the edge states at the 
zigzag edge are at the Dirac energy, shifted away from 
the Fermi level and the superconducting gap. 

Another type of defects which lead to the hybridiza- 
tion of the two valleys are topological defects, like pen- 
tagons or heptagons. These defects can be described at 
long wavelengths as vortices acting on two effective Dirac 
equations with a flux $ = ±(e/i)/(2c), where the two op- 
posite signs ensure time reversal invariance [65l I66j . A 
loop around a defect exchanges the two valleys. A su- 
perconducting gap with opposite signs in the two valleys 
changes sign around such a loop. The resulting system 
can be viewed as two "one fourth" vortices, with each of 
them having a midgap state [571 ISH]) which can be split 
by residual interactions [55]. see also [70]. In the presence 
of a magnetic field, topological defects can pin vortices 
with $ = {eh) /{2c). This configuration leads to an ef- 
fective "half vortex" in one of the two Dirac equations, 
and to a single midgap state, which can be described as 
a Majorana particle. 



Kohn-Luttinger superconductivity and excitonic 
condensation 



The analysis discussed here bears a number of simi- 
larities to the study of excitonic condensation in an ar- 
tificial graphene bilayer when the number of electrons 
in one layer equals the number of electrons in the other 
layer |71H73] , If one assumes isotropic Fermi surfaces in 
both layers, the pairing strength can be described by a 
parameter similar to A,„ in eq. except that the first 
term is negative, and the second term, associated to short 
range interactions, is absent. 

Long range repulsive interactions, which are detrimen- 
tal for the Kohn-Luttinger mechanism of superconduc- 
tivity, are the source of excitonic pairing in an artificial 
bilayer with electrons in one layer and holes in the other. 
An estimate of the effect of long range interactions is 



given in eq. 13 and a similar value was obtained in|73j 
in relation to the excitonic transition. 

An important difference between the Kohn-Luttinger 
mechanism of superconductivity and the excitonic transi- 
tion is the role of a misalignment between the two layers 
and the trigonal distortion of the Fermi surface at high 
carrier concentrations. Superconductivity is due to the 
existence of Cooper pairs where the two carriers involved 
are related by time reversal symmetry, and they are not 
affected by deformations and misalignments of the Fermi 
surfaces (see, however [H]). The electron and the hole 
which are paired and give rise to the excitonic conden- 
sate arise from the same valley, in order for the exciton 
to have zero total momentum. A rotation of the Fermi 
surface, combined with trigonal warping, plays the same 
role as a pair breaking perturbation in a superconduct- 
ing phase. We will not consider here the enhancement of 
the superconducting transition temperature due to fluc- 
tuations of the order parameter, a possibility which has 
been studied in exciton condensation [731 [75]. 



DISCUSSION 

The modulation of purely repulsive electron-electron 
interactions leads to superconductivity with an 
anisotropic order parameter [13]. Pairing by this 
mechanism in isotropic metals is favored when the 
interaction shows a strong momentum dependence, 
V{2kp) 3> V(0) and/or singularities at finite momenta. 
Pairing is be enhanced in anisotropic systems, if the 
interaction which couples regions of the Fermi surface 
with a high density of states is larger than the interac- 
tion at small momenta. Superconductivity induced by 
repulsive interactions which couple different valleys has 
been proposed for the iron pnictides|76l 177] . and it is 
related to the "Van Hove scenario" extensively studied 
in relation to the cuprate superconductors [78] . 

Graphene heterostructures allow for the modification 
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of the long wavelength part of the electron-electron in- 
teraction. The absolute value of the interaction can be 
reduced, making it smaller than the short range repul- 
sion, and the dependence on momentum q can be mod- 
ulated. The lowering of the average value of the long 
wavelength interactions, combined with a significant local 
repulsion, can lead to superconductivity at experimen- 
tally accessible temperatures, for carrier densities above 
5 x 10 13 cm -2 , and interlayer distances of a few nanome- 
ters. The embedding of the device in a dielectric medium 
with a large dielectric constant increases the tendency to- 
wards pairing. 

Superconductivity for an odd momentum pairing state, 
where Ag = — A_g, is unlikely to occur in a single layer 
near a metallic surface 79-82]. The metal can shift the 
chemical potentia of graphene by 0.1-0.4 eV and a charge 
transfer of 10 12 — 10 13 cm -2 is possible. Nevertheless, the 
hybridization with the electrons of the metal will produce 
backscattering, which breaks Cooper pairs and should 
completely suppress superconductivity when the broad- 
ening of the graphene states T is of order T > kBT c . 

In the graphene heterostructures considered here, su- 
perconductivity with a constant gap in each valley, and 
opposite sign in different valleys, is possible. This sym- 
metry implies the existence of rise of triplet Cooper pairs 
with s z = 0. The superconducting phase is gapped, and 
short range scalar disorder is be pair breaking. Surfaces 
and topological defects can induce midgap states. 

We have also found a second superconducting insta- 
bility, induced solely by the intravalley interaction. This 
instability resembles closely the original superconduct- 
ing instability of an isotropic electron gas[T3]. As in that 
case, it occurs at temperatures too low to be accessible 
experimentally. 

Accurate predictions of the value of the critical tem- 
perature of a superconductor are extremely difficult, as 
the value depends exponentially on parameters which are 
not known with great precision. Our estimates depend 
on the electrostatic repulsion between carriers in differ- 
ent layers of the heterostructure, and on the value of the 
short range interaction. The graphene heterostructures 
studied so far show significant effects due to long range 
interactions between layers. The analysis of Coulomb 
drag at high temperatures and high carrier concentra- 
tions seems consistent with models based on the Random 
Phase Approximation, which is also used here. 

In order for superconductivity to occur, the short range 
repulsion must be larger than the long wavelength in- 
teraction. The value of U, the onsite interaction in 
graphene, is difficult to determine experimentally. Cal- 
culations for graphene, aromatic molecules, and conju- 
gated polymers, give values which cluster consistently in 
the range U ^8 — 12 eV, although larger and smaller- 
values have been proposed. For U > 9 eV, supercon- 
ductivity at accessible temperatures can be expected for 
carrier concentrations p > 5 x 10 13 cm~ 2 . Longer range 



interactions, still comparable with the lattice constant, 
will probably enhance the effect of intervalley scattering, 
increasing the value of T c . 

Pairing is very sensitive to the dielectric constant of 
the insulating layers. A lesser source of uncertainty is the 
use of the Fermi energy as the energy scale over which 
pairing takes place. This energy enters in the prefactor 
of the expression for T c . A change in this parameter 
does not alter the balance between long range and short 
range interactions which is responsible for the existence of 
superconductivity. Fluctuations of the superconducting 
phase depress the critical temperature, but their effect 
is not as critical as the details of the electron-electron 
interaction. 
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